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Algebraic and Analytic Aspects of Soliton Type Equations 



Vladimir S. Gerdjikov 



Abstract. This is a review of two of the fundamental tools for analysis of 
soliton equations: i) the algebraic ones based on Kac-Moody algebras, their 
central extensions and their dual algebras which underlie the Hamiltonian 
structures of the NLEE; ii) the construction of the fundamental analytic so- 
lutions (FAS) of the Lax operator and the Riemann-Hilbert problem (RHP) 
which they satisfy. The fact that the inverse scattering problem for the Lax 
operator can be viewed as a RHP gave rise to the dressing Zakharov-Shabat, 
one of the most effective ones for constructing soliton solutions. These two 
methods when combined may allow one to prove rigorously the results ob- 
tained by the abstract algebraic methods. They also allow to derive spectral 
decompositions for non-self-adjoint Lax operators. 



1. Introduction 

We start with three examples of integrable nonlinear evolution equations (NLEE). 



The first one is the well known iV-wave equation ]49j 

(l.i) i[I, Q t ] - i[J, OJ + [[/, Q], [J, Q(x, t)}} = 0, 



L 35 



lim Q(x, t) = 0, 

x — >±oo 



where Q(x,t) is a smooth n x n matrix-valued function, Q(x,t) = —BQ^B and I 
and J arc constant diagonal matrices; B^j — 5ijej, ej = ±1. 



(1.2) 



The second example is the 2-dimensional affine Toda chain [41 
d 2 Q k 



— gQfc+i- 



dxdt 

where we assume that e® n+1 = e® 1 . 



k = 1,. 



The third example belongs to the same family as (1.2) and is of the form: 

n—l 

n, 



.dipk nkd 2 ip k , ■ 'sr* d 



at 



n dx 2 



0=1 



k = 1, 



and k — p is understood modulo n and ipo = ip n = 0. 

The integrability of these equations is based on their Lax representations. This 
means that each of the NLEE (1.1)— (1.3) can be represented as the compatibility 
condition 



(1.4) 



[L(A),M(A)]=0, 
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of two linear matrix differential operators depending on the spectral parameter A. 
Below we will use as Lax operator L(X) 



(1.5) 



L(X)t(}(x, t, A) = ( i— + q(x, t) - A J ) ip(x, t, A) = 0; 



as examples of M(A)-operators we use: 



(1.6a) M{X)*ip = [ i— + V (x,t) + XI ) ip(x,t, A) = Xip(x,t, X)I; 



(1.6b) Mi(A)^ = ( i-jj; + V (x, t) + XV! {x, t) + X 2 V 2 ) ip(x, t, X) = X 2 i>{x, t, X)V 2 a , 



(1.6c) M 2 (X)iP 



i^ + V {x,t) + jV^ 1 (x,t) ) r(r.l.\) 



X 



^(x,t,A)y_ as i; 



where V 2 as = lim x _± 00 V 2 (x, t) and = lim^^ioo V-i(x, t). 

Choosing the form of L(X) in (1.5) we fixed up the gauge by assuming that J 
is constant diagonal matrix and q(x,t) = [J, Q(x, t)], i.e. qjj = 0. 

The system ( fL5| ) with q(x, t) and J 2 x 2-matrices (i.e., q ~ sl(2)) is known as 
the Zakharov-Shabat (ZSs) system; the same system with n x n matrices will be 
referred to below as the generalized Zakharov-Shabat system (GZSs). 

The Lax representation of the N-wave equation is provided by L(X) ( |l.5| ) and 
M(A) ( |l.6a ). If the potentials in these operators are n x n-matrix ones we may 
assume that the Lie algebra underlying the Lax representation is 9 ~ sl(n). The 
set of independent fields Qij(x,t) equals n(n — 1) and may be restricted by the 



involution [49, p2| , 48 
(1.7) 



q{x,t) = Bqt(x,t)B~ 



J = A 



Often by iV-wave equations in the literature people mean eq. ([l]l]) with the invo- 
lution (1.7). Such systems with n = 3 and n = 4 find applications in describing 



wave- wave interactions, sec [48, 49, 35 



The Lax representation of the Z„-NLS eq. (1.3) is provided by (1.5) and (1.6b) 
but with rather specific restrictions imposed on q(x, t) and J: 

n n n 

(1.8) q(x,t)=i^Mx,t)KlS, j = c Y / u k - 1/2 E kk , X = X^W> 



k=l 



k=l 



k=l 



Here and below we will denote by Ejk the nxn-matrices equal to (Ejk) m n — 5j m 5 kn ] 
the indices should be taken modulo n, i.e. E n ^ n+ \ = E n ^ and the constant Co will 
be defined below. 

The affine Toda chain (L2) has several equivalent Lax representations. We 
mention here two of them. Their Lax operators are: 



k=l 



(1.9) L T oda = %~\Y. ^E kk + iXj2 e^-^ 2 E kik+1 , 
and its gauge equivalent: 

(1.10) L Tod . 



k=l 



. d . \ - dQ k 



k=i 



iXKn 
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The corresponding M -operators are of the form (1.6c). Both choices (1.9) and 
( l.lOj ) are not of the form ([PI), but are adjusted to the gr ading o f the Lie algebra 
sl(n,C) we introduce in the next subsection, see formulae ( 1.20| )— (1.25) below. 



The operator LToda (1-10) after a similarity transformation with the constant 



matrix uq, such that u Q 1 K uq — X)fe=i ^E^k can be cast into the form of (1.5) in 
which both q(x,t) and J have a special form: 



(1.11) 



q(x,t) = -i 



^ dx 

fc=i 



Ekk ■ 



where uj = exp(27ri/n). The special form of q(x,t) and J in both (1.8) and ( 1.11 ) 
shows that both models have only n — 1 independent fields. This special form 
can also be made compatible with the structure of the graded and Kac-Moody 
algebras JTT], |||, |l| 



and is best understood with the method of the reduction 



group proposed by Mikhailov 



41 



The idea of the ISM is based on the possibility to linearize the NLEE J53| , [2|, 



^, 13, 49, [|8|, |35 . To this end we consider the solution to the NLEE q(x,t) as a 
potential in L(X) (1.5). In order to solve the direct scattering problem for L(X) we 
introduce the Jost solutions ip±(x, t, A) and the scattering matrix T(A, t) as follows: 

.dip± 



(1.12) 
(1.13) 
(1.14) 



+ (q(x, t) - XJ)ip±(x, t, X) = 0, 



dx 

lim ip±(x, t, X)e 

'. — >±oo 



iXJx 



i. 



T(A,t) = -0+V-(M,A). 

The Jost solutions of L(X) are also eigenfunctions of the operator M (A), 
can use this fact to determine the ^-dependence of the scattering matrix: 

dT 



We 



(1.15) 



dt 



[/(A),T(A,i)]=0, 



which can be easily solved as follows: 



(1.16) 



T(X,t) 



- JfWt 



T(X,0)e 



if(X)t 



By /(A) G \) above we mean the dispersion law of the NLEE; for the N -wave 
system we have /n-w(A) = XI. 

Thus the solution of the NLEE for a given initial condition q(x,t)\ t= o = qo(x) 
can be performed in three steps, see |i~3| : 

1. insert q(x,0) as a potential in L(X) and determine the corresponding scat- 
tering matrix T(A, 0); 

2. Given T(A,0) and the dispersion law /(A) find T(A, t) from eq. ( |l.l6| ); 

3. Given T(X,t) reconstruct the corresponding potential q(x,t) of L(X) which 
will be also the solution of the NLEE. 

Step 2 is trivial. Steps 1 and 3 involve solving the direct and inverse scattering 
problem for (1.5) which can be reduced to linear integral equations. The most 



difficult step 3 provided for the name of the method. The most effective method to 



solve it for operators like L(A) is based on the equivalence to a RHP [ ]45[ ] . 

Along with solving the inverse scattering problem in step 3) we will construct 
also the minimal set of scattering data 1. Indeed the scattering matrix T(X,t) 
has n matrix elements with only one obvious constraint det T(A, t) = 1 while 
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the potential q(x, t) has only n(n — 1) matrix elements. Therefore there must be 
additional interrelations between the matrix elements of T(X,t). 

The analysis of the mapping between q(x, t) and % allows one to interprete it 



as a generalized Fourier transform 



30 




The proof of 



all these facts and the effective solving of the ISP for the GZSs ( |l.5| ) is based on the 
possibility to construct fundamental solutions of (1.5) which are section-analytic 
functions of the spectral parameter A. 



Algebraic structures: Kac-Moody and graded Lie algebras 



Let us now briefly outline the first basic tool inherent in the Lax representation 
- its algebraic structure. Indeed, L(X) and M{X) above are polynomial in A and/or 
1/A whose coefficients take values in some simple Lie algebra g. 

Let us take generic Lax operators in the form: 



(1.17) 
(1.18) 
(1.19) 



i— + U(x,t,X) ) tp(x,t,X) = 0, 



M(X)xjj =\ij t + V(x, t, A) ) ifr(z, t, A) = xl>(x, t, A)^ as (A), 



U(x, t, A) 



k 



U k (x,t)X h 



V(x,t,X) 



J2v k (x,t)X k , 

k 



where the potentials U(x, t, X) and V{x, t, A) are polynomials in A and/or 1/A. Such 
potentials can be viewed as elements of a Kac-Moody algebra gc*. Roughly speaking 
the construction of Qc involves a simple Lie algebra 9 and an automorphism C of 
finite order, i.e. there exist such an integer h that C h = 1. Then we can split g 
into a direct sum of linear subspaces 



(1.20) fl = 

which are eigensubspaces of C, i.e. if 

(1.21) x {k) e Q {k] & 



h-1 



k=0 



C(X {k) )=uj- k X 



(fc) 



where u> = exp(27ri//i). The decomposition ( 1.20 ) satisfies the grading condition: 



(1.22) 



^(k+rn) g „(k+m) ^ 



where the superscript fc + m in g( fc + m ) is understood modulo h. Then the elements 
of the corresponding Kac-Moody algebra Qc have the form: 



(1.23) 



x(x) = xkx{k \ x(k) e Q [k \ 

k<N x 



Obviously due to (1.22) the commutator of any two elements ^(A), Y(X) of the 
form (1.23) will also have the form ( 1.23| ). 

The classification and the theory of the Kac-Moody algebras can be found in 



33, 31 . Their simplest realization can be obtained from a pair (9, C) with a few 



special choices of the automorphism of finite order C, namely: 

a) C = 1; then each of the subspaces ~ g. This leads to a generic GZS 
system if J is real and to a generic CBC system if J is complex. 
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b) C h — 1 where C is the Coxeter automorphism of g and h is the Coxeter num- 
ber. This leads to a CBC system with Z^-reduction and will be used in analyzing 
the NLEE Q and 

c) CV where V is a nontrivial external automorphism of g. Such gradings also 
lead to interesting NLEE but will not be used in this paper. 

The Kac-Moody algebras arc obtained from the constructions a)-c) with addi- 
tional central extensions; they are split into three classes: of height 1 (cases a) and 
b)) and of height 2 and 3 depending on the order of V. 

The potential U(x, t, A) for the TV-wave equations equals [J, Q(x, t)] — XJ belongs 
to a Kac-Moody algebra with g ~ sl(n) and C = 1. The potential U(x, t, A) = 
q(x, t) - A J of the form (Q and fll.llD gives rise to the NLEE (|l]|) and (jt|) is 
related to Kac-Moody algebra of the class b) with g ~ sl(n). The Coxeter number 
then is h = n; the Coxeter automorphism can be realized as inner automorphism 
of the form: 



(1.24) C(X) = C XC„\ C =y> fe £fc/c, u = e 2 ^' n 



k=l 



where C obviously satisfies C n — 1. With this choice of C we can easily check that 
the linear subspaces g' fe ' are spanned by 

(1.25) fl <*> = l.c. {E jJ+k , j,k = 1, . . . ,n} , 

and j + k is considered modulo n. Comparing (L8), ( 1.10| ) with ( 1.25 ) below we 
find that q{x,t) G fl^ and J £ g^K Note that now the condition 6 g^ 

imposes a set of nontrivial constraints on X( k K 

The idea to use finite order automorphisms for the reductions of the NLEE was 
proposed first by Mikhailov [ 41 ] who introduced also the notion of the reduction 



group. The Z„-reduction condition according to |41] is introduced by: 



(1.26) C(U(x,t,Xu)) = U(x,t,X), C(V(x,t, Xu)) = V(x,t,X), 

where we have chosen the simplest possible realization of the 1 n group on the 
complex A-plane: A — > Xuj with u> = exp(2-7ri/n). 

The Kac-Moody algebras, like the semi-simple Lie algebras have an important 
property which ensures the solvability of the inverse scattering problem for L(X) and 
the non-degeneracy of the Hamiltonian structures of the NLEE. While the semi- 
simple Lie algebras possess just one invariant bilinear form (X, Y) = tr (ad x, ad y) 
the Kac-Moody algebras possess a family of invariant bilinear forms: 

(1.27) ((X(X),Y(X)f p) =Rcs n X-p- 1 (X(X),Y(X)), 

for all integer values of p. 

We will need also a central extension of the Kac-Moody algebras g = g c 
where the central element is generated by the 2-cocycle: 

(1.28) u> p (X(x,X),Y(x,X)) = J°° dx(lx(x,X), dY( 2 X) )) iP) - 

This means that each element of g is a pair (X(x, X), cx) where cx is a constant. 
The commutation relation in g is defined by: 

(1.29) [{X{x, A), c x ), (Y(x, A), c Y )} = ([X(x, X),Y(x, X)],u p (X(x, A), Y(x, A))) . 



6 



VLADIMIR S. GERDJIKOV 



Important role for the Hamiltonian formulation of the NLEE is played by the 
dual algebras g*, g* — g* ® c and their splittings into direct sums of Borel-like 
subalgebras. These splittings for g = g + © g_ look like: 



(1.30) 0+ e 

and for the dual g* = g*. 



u k (x)X k 



' -l 

< k= — oo 



u k (x)X k 



(1.31) 



0+ = 1 £ u k {x)\ k 



The co-adjoint orbits of g on g* in fact are isomorphic to the space of coefficients for 
which the NLEE is written. Thus they are natural candidate for the phase space 
of these equations. The freedom provided by the parameter p is directly related to 
the existence of hierarchy of Hamiltonian structures for the NLEE. 



Fundamental analytic solutions 

The second important tool in this scheme is the fundamental analytic solution 
(FAS) of L(X). We will see that using the FAS one is able to: 

- reduce the solving of the ISP for L(X) to an equivalent Riemann-Hilbert 



problem (RHP) for the FAS [45, 48, 



- construct the kernel of the resolvent for L(X) and derive the spectral decom- 
position for L(X) 



26 



18, 30 



- construct the 'squared' solutions of L(X) which all ow t he i nter preta tion of 
the ISM as a generalized Fourier transform (GFT) J2], 



54|, |36|, |25J, |28J, |29 

construct the Green function for the recursion operators A± 



30 



and prove the 

completeness relation for the 'squared' solutions. This property ensures the unique- 
ness of the solution of the ISM §|, |27|, [Tg|, |30 . 



The existence of FAS is ensured by the analytic dependence of both U(x, t, X) 
and V(x, t, X) on A. The properties of FAS depend crucially on the boundary 
conditions imposed on the potential q(x,t). For simplicity here we consider the 
class of potentials q(x,t) that are sufficiently smooth functions of x and tend to 
zero fast enough for x — > ±00 for any fixed value of t. 

The FAS for the Zakharov-Shabat system (i.e. g ~ sl(2)) can easily be con- 
structed due to the fact that each of the columns of the Jost solutions 



(1.32) 



L(A)V>±0M,A) =0, 



lim e 

x — *±oo 



iJXx 



ip±(x,t,X) = 1, 



allow analytic extension either for A G C+ or for A £ C_, sec Q{. 

If we analyze the analyticity properties of the Jost solutions tp±{x, t, X) related 
to algebras of higher rank one finds that only the first and the last columns of 
ip±(x,t,\) allow analytic extensions off the real A-axis. An important result of 
Zakharov and Manakov [49, 48] consisted in showing that a FAS for the GZS with 



g ~ sl(n) and real- valued J can be constructed by taking proper linear combinations 
of the columns of ip±(x, t, A). 

The construction is more complicated for the Caudrey-Beals-Coifman (CBC) 
systems when the eigenvalues of J are complex S E$L ||| . The generalization of this 
construction for CBC systems related to any simple Lie algebra g was done in |30|. 
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We make attempt to outline the construction and the properties of each of 
these tools. Then we show how the FAS can be used to construct the kernel of 
the resolvent of L(X) and to exhibit its spectral properties and the structure of its 
discrete spectrum. Finally we illustrate how these tools can be used in the analysis 
of the NLEE and their fundamental properties and finish with some conclusions. 

Both these aspects are rather broad; they have been widely discussed in hun- 
dreds of papers. Therefore inevitably the list of references consists mainly of reviews 
and monographs and bears an illustrative character. The thorough reader is advised 
to consult also the papers referred to in these references. 



2. Construction of the FAS 

Preliminaries: Jost solutions and scattering matrix 



The direct and the inverse scattering problem for the Lax operator (1.5) will 
be done for fixed t and in most of the corresponding formulae t will be omitted. 

The crucial fact that determines the spectral properties of the operator L is 
the choice of the class of functions where from we shall choose the potential q(x). 
Below for simplicity we assume that the potential q{x) is defined on the whole axis 
and satisfies the following conditions: 

C.l: By q(x) € 9Jts we mean that q(x) possesses smooth derivatives of all 
orders and falls off to zero for \x\ — > oo faster than any power of x: 

lim \x\ k q{x) = 0, Vfc = 0,1,2,... 

x — >±oo 

C.2: q{x) is such that the corresponding operator L has only a finite number 

of simple discrete eigenvalues. 
Below we will use along with L^p(x, A) = also the following equivalent formu- 
lations of the system fll.E 



(2.1) i^L + q(x,t)Z(x,\)-*[J,fa,X)]=0, Z(x,\) = iP(x > \)e iXJx , 

(2.2) i^-tp(x,X)q(x,t) +\i>(x,\)J = 0, ${x,X) = (ip(x,X))-\ 

(2.3) i^--i(x,X)q(x,t)+X[i(x,X),J]=0, £(x, A) = e' lXJx ^(x, A). 

ax 

where by 'hat' we denote the inverse matrix, X = X^ 1 . The Jost solutions £±(x, A), 



X±(x,X) and £±(x, X) for the systems ( J2 . l| ) — ( 2.3) can be introduced by: 

lim £±(.t, A) = 1, lim if}±(x, X)e~ iXJx = 1, lim £±(x,X) = l, 

X— >±oo X~ >±00 >±OC 



in analogy to (1.13); then their scattering matrices are: 

T 2 (X)=e iXJx (^(x,X))-H+(x,X)e- iXJx =T(X), 

T 3 (A) = 1 jj + (x 1 X)(^(x,X))- 1 =T- 1 (A), 

T 4 (A) = e iXJ *i + (x, X)(t(x, X))~ l e- iXJx = T~\X), 

Below we will consider two specific reductions of the Lax operator: the GZSs 
with Z2-reduction: 

(2.4) B{U i <(x,t 1 eX*))B- 1 = U(x,t,X), B 2 = 1, e = ±1. 
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The first possible choice for B = diag(ei, . 
classical iV-wave equations |40| . with 



. , e„), €j = ±1 with e = 1 leads to the 



(2.5) 



Qkj( x ^) = ^jQjk(x,t), J = diag(ai, . . .,a n ), 



a k = ea k . 



Since all eigenvalues of J are real (e = 1), or purely imaginary (e = —1), the Lax 
operator becomes a GZSs. The second choice for B: 



(2.6) 



B 



i 



fc=i 



will be used in combination with the Z„-reduction: 
(2.7) CoO/M^A))^ 1 = U(x,t,X), 



1, 



C" = 1. 



which leads to the CBC system. For the sake of convenience in doing the spectral 
problem of CBCs we choose Co = J2k=-\ ffis.fc+ jj th en L(X) has the form ( |l .5| ) with 
diagonal complex-valued J given by (1.8) or (1.11) where Co = 1 (resp. Co = i) if 
e = 1 (resp. e = —1). Both Lax operators will have similar spectral properties. 



In solving the NLEE (1.2) and fll-q ) we will need to apply both reductions 
( p~i| ) and (2.7) simultaneously. An attempt for classification of the Z 2 -reductions 
is made in [ |23[ . 



The FAS of the GZSs with Z 2 -reduction. 

Let us outline without proofs the construction of the FAS for the GZSs with 
real J, see [ 49 , |48| , ||, ^| . For dcfiniteness we assume that the real eigenvalues 
of J are pair- wise different and ordered as follows: 



(2.8) 



J = diag (ai, . . . , a n ), ai > a 2 > 



Proposition 2.1. Let the potential of (1.5) q(x) £ Wis satisfies conditions 



(C.l), (C.2) and (2.1). Then: 



a) the Jost solutions £±(x,A) and $,±(x,A) of (2.1), (2.i) exist and are well 
defined functions for A £ M; 

b) the matrix elements of the scattering matrix T(X) and its inverse T(A) are 
Schwartz-type functions of A for A G M. 



Remark 2.2. The proposition 2.1 concerns the Jost solutions as fundamental 
solutions. One can prove that the first and the last columns f^± (x, A) and £± ' (x, A) 
of the Jost solutions allow analytic extension with respect to A as follows: 

Column Z l + ] (x,\) £, l + ] {x,\) £™(ar,A) C- ] (x,X) 

Analytic for A e C_ A G C+ A e C+ A G C_ 

Analogously the first and the last rows of the Jost solutions (x, A) and £± (x, A) 
allow analytic extension with respect to A as follows: 



Row 

Analytic for 



e£ ] (s,A) i [ : ] (x,X) B ] (x,X) t ] (x,X) 



A G 



A G C_ 



A G C_ 



All the other columns of £±(x, A) and rows of £±(x, A) are defined only for A G 
and do not allow analytic extensions off the real axis. 
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We start by explaining the construction of the FAS x ( x >^) or rather of the 
solutions 



(2.9) 



£ ± (x,\)=x ± (x,\)e 



iXJx 



45 



to eq uat ion (2.1) which allow analytic extensions for A G C±. Skipping the details 



, 48, 49 1) we formulate the answer and determine £ + (x, A) as the solution 



of the following set of integral equations: 

(2.10a)£±(iz, A) = 5 l3 + i f dye-^-^*^ £ q ip (y)^(y, A), i > j; 

J — OO „_i 



P =i 



(2.iob)c+(x,A) = i / dye-^~ a ^-y^^ P (yK 3 (y^l i<r, 
J °° p=l 

Analogously we define £,~ (x, A) as the solution of the set of integral equations: 

(2.11a)^(x,A)=i I* dye-^-^-rt^M&'V' l> X 
Jo ° P =i 

(2.Ub)%(x,X) =5 ij +i [ X dye-^-^-^^M&'V' f ^'i 



Theorem 2.3. Let q(x) G Ms satisfies condit ions (C.l), (C.2) and let J sat- 
i sfy (!Ti\ ). Then the solution £ + (x, A) of the eqs. ( 2.1(\) (resp. (;~(x,X) of the eqs. 
( 2.11 )) exists and allows analytic extension for A G C+ (resp. for A G C-). 



Remark 2.4. Due to the fact that in eq. ( 2.10 ) we have both oo and -co as 
lower limits the equations are rather of Fredholm than of Volterra type. Therefore 
we have to consider also the Fredholm alternative, i.e. there may exist finite number 
of values of A = A^ G C± for which the solutions ^(x, A) have zeroes and pole 
singularities in A. The points A^ in fact are the discrete eigenvalues of L(X) in C±. 



The reduction condition ( |2.4| ) with e = 1 means that the FAS and the scattering 
matrix T(A) satisfy: 

(2.12)5 ( X +(x, A*)) f B- 1 = ( X -(x, A))" 1 , B (T(A*)) f B^ = (T(A))- 1 . 

Each fundamental solution of the Lax operator is uniquely determined by its 
asymptotic for i -> oo or i -> — oo. Therefore in order to determine the linear 
relations between the FAS and the Jost solutions for A G K we need to calculate 
the asymptotics of FA S for x — ► ±00. Taking the limits in the right hand sides of 
the integral equations ( 2. 1C ) and (2.11) we get: 



(2.13a) lim &(x, A) = 5u, for i > j; 

X^ — OG J 

(2.13b) Jim _j7.( x ,\)=S ij , for i < j; 



X— — OO 



lim &(x, A) = 0, for i < j; 
lim £7. (x, A) = 0, for i > j; 

-r.- — >no J 



This can be written in compact form using (2.9): 

(2.14) X ± (x, A) = j>-(x, X)S ± (X) = i/,+ (x, X)TT(X)D ± (X), 
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where the matrices S ± (X), D ± (X) and T ± (X) are of the form: 



(2.15a) S+(A) 



(2.15c) S~(X) 



/1S+... 
1 ... 


Sin \ 

°2n 


\0 ... 




diag (D+,D 


+ 

2 , . . . 


/ 1 . 

s 2 - x 1 . 


.. o\ 

.. 


V ^nl ^712 • 


.. 1/ 



T+(A) = 



1 T+ 

1 J 12 
1 



2ft \ 
T+ 1 



\0 ... 1 / 

IT (A) =diag(Z>r, J D 2 - 
/ 1 ... o\ 



,D-), 



T-(X) 



T 2 'i 1 



1/ 



Let us now relate the factors T ± (A), 5 ± (A) and I? ± (A) to the scattering matrix 
T(A). Comparing ( |2.14f ) with fll.l^) we find 

(2.16) T(A) = T-(A)Z?+(A)S ,+ (A) = r + (A)D-(A)5"(A), 

i.e. T ± (A), S' ± (A) and D ± (A) are the factors in the Gauss decomposition of T(A). 

It is well known how given T(A) one can construct explicitly its Gauss decom- 
position, see the Appendix |a]. Here we need only the expressions for D ± (X): 

m ^ X) D -(X) = m ""i+ l(A) 



(2.17) 



m+(xy 



m -(A) 



where rrv^ are the principal upper and lower minors of T(A) of order j. 

Corollary 2.5. The upper (resp. lower) principal minors rrif (A) (resp. mJ(X) 
ofT(X) are analytic functions of X for X G C+ (resp. for X G C-). 



Proof. Follows directly from theorem |2.3| , from the limits 
(2.18) 



lim £+(x,A) =D+(A), 

x — *oo JJ J 



lim LAx,X) = Dj (A), 

X — > OO J J J 



and from (|2.15b|) and Q2.17|) . 



□ 



Corollary 2.6. The following relations hold: 



a) lim £ (x, A) 

A — >oo 



6J lim mf (A) = 1. 



Proof, a) follows from the integral equations ( 2.10| ), (2.11) taking into account 
t hat t he int egrand s in their right hand sides vanish for A — > oo. b) follows from a), 
( fjf ) and ( |2.15b| ). □ 

In what follows we will also assume that the set of minors (A) have only 
finite number of simple zeroes located at the points 

(2.19) 3 = {A±GC±, ./ 1 V.} 

Generically each of the A^ can be a zero of a string of minors, e.g.: 

(2-20) m+(A) = (A - A+)m+ 3 + O ((A - A+) 2 ) , 

for 1 < Ij < Fj < n. Let us introduce the quantities bjk as follows: 

1 if A+ is a zero of to^(A); 
if A^~ is not a zero of m^"(A). 



(2.21) 
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and note that the reduction (2.4) means that the Gauss factors of T(A) satisfy 
( £ = 1): 



(2.22a) 
(2.22b) 



B 



(£ + (A*)) B- 1 = S~(X), B (f + (X*j) B~ l =r~(A), 



D+(X*)) = D~(X). 



The relations (2.22a) are strictly valid only for A £ K while (2.22b) together with 
( [2.15b ) and ( 2.17 ) leads to the following constraints on the minors to^(A): 

(2.23) (m+(A*))*=m-_ fc (A). 

Thus if A^T is a zero of m^(X) then A^ = (A^T )* is a zero of m~_ k (X). 

The FAS of the CBCs with Z„-reduction. 

The crucial dif fere nce with the Z2-case treated above consists in the fact that 
now J is given by (1.8) or (1.11) and has complex eigenvalues. Skipping the details 
(see JB], ^]) we just outline the procedure of constructing the FAS. 

First we have to determine the regions of analyticity. For potentials q(x) sat- 
isfying the conditions (C.l) and (C.2) and subject to the Z n -reduction (2/7) these 
regions are the 2n sectors fi„ separated by the rays l v o n wh ich Im X(aj — a^) = 0. 
We remind that if we assume also the Z2-reduction (2.£) with Cq = eco then 
ak — cqlj^ 1 / 2 . Then the rays l w are given by: 

it(y - 1) 



(2.24) 



l u : arg(A) = <j) + 



,2n, 



0. Thus the 



where <Pq — 7r/(2n) only if e = 1 and n is odd; in all other cases <f>Q 
neighboring rays l v and close angles equal to ir/n. 

The next step is to construct the set of integral equations analogous to (2.10) 
whose solution will be analytic in To this end we associate with each sector Q, v 
the relations (orderings) > and < by: 



(2.25) 



i > j 

V 

i < j 



if 



Im \(cii 
Im A(ai 



aj) < for A e fl v , 
aj) > for A e Cl v . 



Then the solution of the system ( p. 10 ) 



(2.26a)$(z,A) = *y+i 



(2.26b)&(x,X)=i 



[ X dye-^- a ^ x ^ ]T qM&iv, A), 

J OO i 



dye-iH^-a^-y) q ip (y)Q(y, A), i > j; 
P =i 

h 

i<j; 

P =i 

will be the FAS of the CBCs in the sector f2„. The asymptotics of £"(a;, A) and 
A) along the ray l v can be written in the form: 

(2.27a) lim e iXJx ^(x, Xe M ) e - lXJx = 5+ (A), A € l„, 

> - — OO 

AXJx 1 /„ \ ^—iO\^—i\Jx 



(2.27b) 
(2.27c) 
(2.27d) 



lim e" 

X — > — oo 



lim 



*i v ~ l {x,Xe Y 
c r _i (a;,Ae- 40 )e- 



lXJx C{x,Xe M )e- lXJx = T~D+(X), 



SuW, 



lim e 

x — *oo 

i\jx£v— 1 \ „ — iO \ „— iX Jx D~ (A) 



A e l u , 
X G l„, 
X e l v , 
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where the matrices S+, Tj~ (resp. S~ , T~) are upper-triangular (resp. lower- 
triangular) with respect to the ^-ordering. As in the previous case they provide the 
Gauss decomposition of the scattering matrix with respect to the ^-ordering, i.e.: 



(2.28) 



T„(A) = T-(A)I>+(A)S+(A) = T+(X)D-(X)S-(X), 



A e l v . 



More careful analysis shows [ pOf that in fact T V (X) belongs to a subgroup (3„ of 
SL(n,C). Indeed, with each ray l v one can relate a subalgebra q„ C sl(n,C). 

If Z„-symmetry is present each of these subalgebras g v is a direct sum of sl(2)- 
subalgebras. Each such s/(2)-subalgebra can be specified by a pair of indices (k, s) 
and is generated by: 



(2.29) 



kk 



Ei. 



f 



(k,s) 



k< s. 



Then the scattering matrix T U (X) will be a product of mutually commuting matrices 

T!) k ' s) of the form: 

(2.30) 

= 1 + (a±, JX) - l)E kk + (a-.,AX) - 1)E SS - bZu 3 (X)E ks 



where k< s, with only 4 non-trivial matrix elements, just like the ZS (or AKNS) 

V 

system. 

The Z„-symmetry imposes the following constraints on the FAS and on the 
scattering matrix and its factors: 



r- 2 (*,A), 



C TAXu)C 



T„_ a (A), 



(2.31a) CoffaXu)^ 1 

(2.31b) CoStiX^Co 1 = fiSL a (A), C^Au;)^- 1 = £>*_ 2 (A), 

where the index ^ — 2 should be taken modulo 2n. Consequently we can view as 
independent only the d ata o n two of the rays, e.g. on l\ and l2 n = lo', ah the rest 
will be recovered from ( 2.31 ) . 

If in addition we impose the Z2-symmetry ( |2.4| ), (2.6) with e = — 1 then we will 
have also a k = ito^ 1 / 2 and: 

(2.32) 

sr^-A'ijtr 1 = (e +1 -"(x,x))-\ Bistix^B- 1 = (S^ +1 _ u (X))-\ 

and analogous relations for Tj^(X) and D^(X). Another interesting subcase takes 
place for even values of n and ^-reduction (2.4), (2.6) with e = 1; then a k — oj k ~ 1 / 2 
and the FAS satisfy: 

(2.33) 

5(0, A*))^- 1 = (e n+1 -»(x,X))-\ BiStiX^B- 1 = (S? n+1 _ v (X))-\ 



In both cases the rays l v are defined by ( 2.24 ) with (f>o — 0. The pairs of indices 
{k u ,m v } specifying the imbeddings of the sZ(2)-subalgebras related to the ray l v 
are defined as follows: 



(2.34) 



a) for e 

b) for e 



1 



k v + m v = 
k,, + m„ 



f 2 — v{ mod n), 
i/( mod n), 



One can prove also that D+(X) (resp. D u (A)) allows analytic extension for 
A G il„ (resp. for A G compare with corollary Another important fact is 

|3§ that D+(X) = D~ +1 (X) for all A G Cl v . 
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The inverse scattering problem and the Riemann-Hilbert problem. 

The next important step is the possibility to reduce the solution of the ISP for 
the GZSs to a (local) RHP. Indeed the relation ( |2.14[ ) can be rewritten as: 

(2.35a) £+(x,t,X) = £-(x,t,X)G(x,t,X), AeR. 

(2.35b) G{x,t,X) = e -i(A^-/(A)t) G < o ( A ) e i(AJx-/(A)t) ! 

(2.35c) G (A) = S-(X)S + (X) ; 

in other words the sewing function G{x, t, X) satisfies the equations: 

(2.36) i^-X[J,G(x,t,X)}=0, i^ + [f(X),G(x,t,X)}=0, 
dx at 

Here /(A) G f) determines the dispersion law of the NLEE. Together with 



(2.37) 



lim t*(x,\) = 1, 



eq. ( 2.3 5| ) is known as the RHP with canonical normalization. 

Theorem 2.7 (|45|). Let £ + (x,t, A) and £~(x,t,X) be solutions to the RHP 



(2.35), (2.31) allowing analytic extension in A for A G C± respectively. Then 
X ± {x, t, A) = ^(x, t, X)e lXJx are fundamental analytic solutions of both operators 
L and M , i.e. satisfy eqs. (1-6) with 

(2.38) 



q(x,t)= lim \ (j - £ ± (x,t,\)ji ± {x,t,\)) 



Proof. Let us assume that ^(x,t, A) are regular solutions to the RHP and 
let us introduce the function: 



(2.39) 



g ± (x, 1, A) = i-^r^fa *. A ) + ^(x, t, X)J£ ± (x, t, A). 



If ^(x, t, A) are regular then neither ^(x,t, A) nor their inverse ^(x,t, A) have 
singularities in their regions of analyticity A G C±. Then the functions g^(x,t, A) 
also will be regular for all A G C±. Besides: 

(2.40) 

The crucial step in the proof of 



lim g (x,t,X) = lim g (x,t,X) = XJ. 

A — *oo A — >oo 



52 



is based on the chain of relations: 



9 (x,t,X) = i 



(2.35) .d(£ G) 



dx 



Gf(i, t, A) + XCGJGC(x, t, A) 



.dG 



i^-i-(x,t,\)+£- ( i^G + XGJG{x,t,X) ) £-(*,*, A) 



dx 



dx 



^ i^-C(x,t,X)+C (x[J,G]G + XGJG(x,t,X)) i~(x,t,\) 
dx V / 

= i^rC(x,t,\) + \CJi-(x,t,\) 

(2.41) = g-(x,t,X), AgK. 

Thus we conclude that g + (x,t, A) = g~(x,t, A) is a function analytic in the whole 
complex A-plane except in the vicinity of A — > oo where g + (x,t, A) tends to A J, 
( 2.40| ). Next from Liouville theorem we conclude that the difference g + (x, t, A) — A J 
is a constant with respect to A; if we denote this 'constant' by —q(x,t) we get: 

(2.42) g + (x,t,X)-XJ = -q(x,t). 
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It remains to remember the definition of g + (x,t, A) (2.39) to find that £ (x, t, A) 
satisfy (2.f), i.e. that x ± (x,t, A) is a fundamental solution to L. The relation 
between q{x 7 t) and ^(x^A) ( 2.38| ) can be obtained by taking the limit of the 
left-hand sides of (2.42) for A — * oo. 

Arguments along the same line applied to the functions h (x,t, A) 



(2.43) 



/i±(x, i, A) = i^f^s, t, A) - ^(x, f, X)fW^(x, 1, A), 



can be used to prove that x^( x : ^ A) are fundamental solutions also of the operator 
M; equivalently it satisfies (V'(x,t, A) = V(x,t, A) — /(A)): 

(2.44) ,^_ + y'( a;; i ; AK ± ( a :,t,A) + [/(A)^ ± ( a :,t,A)] =0, 

and one finds that h + {x,t, A) = h~{x,t, A) is a function analytic everywhere in C 
except at A — > oo where it has a polynomial behavior of order N — 1. Denoting the 
polynomial as V(x, t, A) we derive ( 2.43| ). 

To conclude the proof of the theorem we have t o acc ount for possible zeroes 
and pole singularities of ^(x,t, A) at the points 3 ( |2.19 ). Below we derive the 
structure of these singularities which is such that they do not influence the functions 
g ± {x 1 1, A) and h ± (x, t, A). The theorem is proved. □ 

The analyticity properties of m k (A) allow one to reconstruct them from the 
sewing function G(A) ( 2.35c ) and from the locations of their zeroes through (see 
Appendix H) : 



(2.45) S) fc (A) = 

where 
(2.46) 



dp 



2tt« J-oo M - A 



In 



1, 2, 
1, 2, 



iV 



v - A - A+ 
5> fc ln 



A- A- 



S fe (A) 



_/lnm+(A), 

lnm „-fe( A ) 



A e 



A G 



One can view 55fc(A) as generating functionals of the conserved quantities for 
the related iV-wave-type equations; the relevant expressions for them in terms of 
the scattering data can be obtained from the right hand sides of ( 2.45 ). 

Quite analogously we can treat also the CBCs with Z„-symmetry. More pre- 
cisely, we have: 



(2.47) 



A e l u , 



G v (x,t,X) = e- iXJx+if ( X)t G 0}U (X)e iXJx - ifmt , 



The collection of all relations (2.47) for v = 1,2, 



Go,, (A) = S-(X)S+(X) 
, 2n together with 



t=o 



(2.48) 



lim C{x,t,\) 



can be viewed as a local RHP posed on the collection of rays £ = {l u }'^L l with 
canonical normalization. Rather straightforwardly we can reformulate the results 
for the GZSs for the CBCs and prove that if (x, A) is a solution of the RHP ( gig ) , 
( gig ) then x"{x, A) = $ v (x, X)e lXJx satisfy the CBC with potential 



(2.49) 



q{x,t)= lim ( J -Z v (x,t,X)j£' / {x,t,X)) 
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We finish this subsection by formulating the dispersion relations for the functions 
In tyi v j. (A) which allows us to reconstruct them from their analyticity properties: 



(2.50) 
lnTO ifc( A ) 



In 

E 

ij=i 



2ni 



dfi 
/i — A 



hi 



k 



N 



EE^ ln 



A- A 



n =lj=l 



where A G and the superscript r\ in the integrand shows that we should use the 
ordering characteristic for the sector tt n ; b v kj a re the analogs for bkj (2.21) in £l v . 

Both dispersion relations ( p. 45 ) and ( 2.50| ) can be used to derive the so-called 
trace identities (see [48, 13|) for the GZSs and CBCs respectively. Indeed, S)fe(A) 
and lnm^ fe (A) allow asymptotic expansions 



(2.51) 



S) fc (A) = 5> 



lnm+ fe (A) 



E M S A " 



The expansion coefficients 2)^ and Af^ / are local integrals of motion, i.e. their 
densities depend only on q(x, i) and its derivatives with respect to x. Their explicit 
calculation is done via recurrent procedure. We illustrate it by the two first integrals 
of motion of the Z„-NLS equation (1.3): 



-(1) 


= -/ 


1,1 


2lu J_ 


fl 






= 2^] 



■0p 



dx 



- E 

3 f-f 

p+fc+Z— n 



tp p tpkipi{x,t) 



One can also expand the right hand sides of the dispersion relations ( |2.45| ) and 
( |2.50|) over the inverse powers of A which allows to express 35 k and Mjfl also in 
terms of the scattering data of GZSs and CBCs. 



The dressing Zakharov-Shabat method 

One of the most fruitful ideas for the explicit constructing of the NLEE's solu- 
tions is based on the possibility starting from a given regular solutions £ (x, t, A) 
to the RHP to construct new singular solutions £ (x,t, A) having zeroes and pole 
singularities at the prescribed points A^ G C± . The structure of these singularities 
are determined by the dressing factor Uj(x,t, A): 

(2.54) ^(x,t,X) = u j (x,t,X)£(x,t,X) U r 1 _(*), 

which for the SX(n)-group has a simple fraction-linear dependence on A: 

(2.55) Uj(x,t,X) = t + (c j (X)-l)P j {x,t), 

(2.56) uj]_ = lim uAx,t, A). 



9(A) 



A -A- 
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Here Pj (x, i) is a projector Pf = Pj which for simplicity is chosen to be of rank 1; 



then it can be written down as: 



\nj){mj\ 



(nijlnj) ' 

where the bra- and ket- eigenvectors (m,j\ and \rij) are the 'left' and 'right' eigen- 



(2.57) Pj(x) 

where the bra- and ket- 
vectors of the projector. 

From (2.54) there follows that the dressing factor u(x,t, A) satisfies the equa- 
tion: 

du 

(2.58) i— — h q(x, t)u(x, t, A) — u(x, t, X)qo(x, t) — A[J, u(x, t, A)] = 0. 

ax 

The main advantage of the dressing method is in the fact that one can determine 
the x and i-dependence of (rrij\ and \rij) through the regular solution Xo ( x > 1> A) as 
follows: 

(2.59) K) = x+.(x,t)\r$), K| = (m^lxoj^t), Xofat) = Xo(x,t,\f) 
or equivalently these vectors are solutions to the equations: 

(2.60) i-M + ^°)(a: > t,At)|» J .) = J i^- + V^(x,t,Xf)\n j ) = 0, 

ax at 

(2.61) ifei - ( mj \uW(x,t,\j) = 0, »-M - ( mj \V^(x,t,XJ) = 0, 

(2.62) U(°\x,t,X)=q (x,t)-\J, V^(x, t, A) = V(x, t, X)\ q=qo . 

Here qo(x,t) is the potential corresponding to the regular solut io ns Xp (a;, t, A) to 
the RHP and V^(x, t, A) is obtained from V(a;, i, A) (see ( 3.35 ), ( 3.36 )) replacing 
g(x, t) by go(£, This construction is well defined also in the case when Xo ( x i ^) 
are singular solutions to the RHP, provided they are regular for A = A^ . 

If q(x,t) is the potential corresponding to the singular solution x \ x , ^ A) then: 

q(x,t) = qo(x,t) + lim A(J — Uj(x,t, X)Juj(x,t, A)) 

A — >oo 

(2.63) = q (x,t) - (A+ - \J)[J, Pj(x, t)]. 

Thus starting from a given regular solution of the RHP (and related solution 
qo(x,t) to the NLEE) we can construct a singular solution to the RHP and a new 
solution q[x, t) of the NLEE depending on the Xj 1 and on the eigenvectors of Pj(x). 
If we start from the trivial solution qo(x,t) = of the NLEE then we will get the 
one-soliton solution of the NLEE. Repeating the procedure N times we can get the 
iV-soliton solution of the NLEE. 

With the explicit formulae for Pj(x) and using ( 2.54 ) we can establish the 
relationship between the scattering data of the regular RHP and the corresponding 
singular one. The dressing factor Uj(x, A) is determined by the constant vectors 
(m,j\ and \n®) can not be quite arbitrary. The condition that q(x) vanishes for 
x — > ±oo requires that if (n^) 3 = for all 1 < s < Ij and Fj < s < n then also 
(mj) s — for all 1 < s < Ii and F\ < s < n. Thus we derive that: 

(2.64) lim Pj{x) =E IjZj , lim Pj (x) = E FjF . , 

x — *oo x — y—oo 

and therefore 

(2.65) u j , + (X)=l + (c j (X)-l)E IjIj , Uj ,-(X) = l + ( Cj (X)-l)E FjFr 
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The interrelation between the Gauss factors of the corresponding scattering matri- 
ces are: 



(2.66) S ± (X) = u jt _(A)Sf(A)t£i(A), T±(A) = % -,+(A)T ± (A)^(A), 
and 

(2.67) D ± (X) = u A+ (A)ltf (A)i£l(A). 



Comparing these last relations with ( 2.17 ) we find for the principal minors of T(A) 
and T (A): 



(2.68a) m+(A) 
(2.68b) mj(X) 



X-X] 
X-X- 
A-A, 
A-Aj 



-m£ s (X), for Ij <s<Fj, X G C+ U : 



■mQ S (A), for n 



Fj < s < n — Ij, 



and ra^(A) = m t s (A) for the other values of s. Thus the result of the dressing is 
that the string of upper principle minors mj(X), Ij < s < Fj acquire simple zero 
at A = A^~ while the string of lower principle minors mj(X), n — Fj < s < n — Ij 
acquire simple zero at A = A^ . 

Obviously if we impose on L(X) the ^-reduction then we should restrict also 
the dressing factor by: 



(2.69) 



B (u(x, t, eA*)) f B- 1 = u(x, t, A). 



The ansatz (2.55) satisfies ( ^.69 ) if A^ = e(A^)* and the vectors \noj), (^Oj'l are 
related by: 

(2.70) (m 0j HS| ? 4). 

If we impose the Z„-reduction ( ^.7] ) then u(x, i, A) must satisfy: 

(2.71) Couix^uXjCo 1 = u(x,t,X). 

Such conditions require generalizations of the ansatz ( |2.55 ) [41 



(2.72) 



,{x, t, X) = 1 + ]T ( Cj (cj s X) - 1) C*P 3 {x)C 



A slightly different approach treating also multi-soliton solutions of the Z„-symmetric 
NLEE is given in ||. 

Up to now we dealt with the algebra g ~ sl(n, C). Treating the other simple Lie 
algebras (orthogonal or symplectic) needs additional care especially in constructing 
the dressing factors 51, 23]. 

In fact Uj(x , A) ( ^.55 ) must be an element of the corresponding group. From 
the ansatz (2.55) it follows that Uj(x,X) belongs to GL(n,C), but one can always 
multiply u(x, X) by an appropriate x- and i-independent scalar and to adjust its 
determinant to 1. Such a multiplication goes through the whole scheme outlined 
above but is adequate only for the sl(n, C) case. However the ansatz ( [2.55 ) can not 
be used, e.g. for the case so(n,C). The adequate ansatz is formulated below 
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Theorem 2.8. Let g ~ B r or D r and let the dressing factor u(x, A) be of the 
form: 

(2.73) 

Uj (x, X) = i+ (c,(A) - l)Pj(x) + (c-\\) - l)P-j(x), P H = SvPjSv 1 , 

where So is introduced in (A. lip and Pj(x) is a rank 1 projector (2. 51). Let the 
constant vectors \hq) and (mo| satisfy the condition 

(2.74) (m \S\m Q ) = (n \S\n ) = 0. 

Then uj(x, A) ( 2. 7c\ ) satisfies the equation ( 2.5^) with a potential 



(2.75) 



q(x) = q (x) - (A+ - A - )[J,pj(x)}, pj(x) = Pj(x) - P-j(x). 



PROOF. Due to the fact th at Xo ( x > ^) take values in the corresponding orthog- 
onal group we find that from (2.74) it follows (m\S\m) — 0, (m\JS\m) = and 
analogous relations for the vector \n). As a result we get that 

(2.76) Pj (x)P. j (x) = P.j (x)P 1 (x) = 0, P 3 (x)JP-j (x) = P_j (xpPj (x) = 0. 



Let us now insert ( 2.73 ) into ( 2.58| ) and take the limit of the r.h.s idc of ( 2.58 ) 
for A — » oo. This immediately gives eq. ( |2.75|) . In order that Eq. ( ^.58| ) be satisfied 
identically with respect to A we have to put to also the residues of its r.h.side 
at A -> A+ and A -> Aj. This gives us the following system of equation for the 
projectors Pj(x) and P-j(x): 



(2.77) 
(2.78) 



.dP, 



+ q W P iW ~ P ^ x ^(x) - \J{J,Pj(x)} = 0, 



i-^f- + q{x)P- j (x) - P- 3 (x)q (x) - At [J, P^{x)] = 0, 



where we have to keep in mind that q is given by (|2.75| ). Taking into account (|2.76[) 
and the relation between Pj(x) and P-j(x) eq. ( 2.77|) reduces to: 

(2.79) 



dP 

i-£ + ko(x), Pj(x)\ + XjPj(x)J - \+JP 3 {x) - (\j 



X+)P j (x)JP j (x)=0. 



One can check by a direct calculation that ( [2.57 ) satisfies identically (2.79). The 
theorem is proved. 

□ 



3. The resolvent and spectral properties of GZSs and CBCs 

The FAS x ± ( a; ! ^) °f L(X) allows one to construct the resolvent of the operator L 
and then to investigate its spectral properties. By resolvent of L(X) we understand 
the integral operator R(X) with kernel R(x,y, A) which satisfies 

(3.1) L(X)(R(X)f)(x) = f(x), 

where f(x) is an n-component vector function in C™ with bounded norm, i.e. 
C 00 dy(f T {y)f( V ))<^. 

From the general theory of linear operators 12, 46 1 we know that the point 



A in the complex A-plane is a regular point if R{X) is a bounded integral operator. 
In each connected subset of regular points R(X) is analytic in A. 
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The points A which are not regular constitute the spectrum of L(X). Roughly 
speaking the spectrum of L(X) consist of two types of points: 

• i) the continuous spectrum of L(X) consists of all points A for which R(X) is 
an unbounded integral operator; 

• ii) the discrete spectrum of L(X) consists of all points A for which R(X) 
develops pole singularities. 

Let us now show how the resolvent R(X) can be expressed through the FAS of 
L(X). Indeed, if we write down i?(A) in the form: 

/>oo 

(3.2) R(X)f(x) = / R(x,y,X)f(y), 



the kernel R(x, y, A) of the resolvent is given by: 

(3.3) R( X ,y,X) = (* + ^ y ^\f rX ' G ^ + > 
y 1 K ,y ' ' \ R~(x,y, A) for A G C~, 

where 

(3.4) R ± (x,y,X) = ±i X ± (x,X)e ± (x~y)x ± (y,X), 

kg 

e ± (z) = 8( T z)ii Q - 9{±z){t - no), n = £ e ss , 

s=l 

where ko is the number of positive eigenvalues of J; namely: 

(3.5) oi > a 2 > • • • > a ko > > a ko+1 > ■■■ > a n . 

Due to the condition tr J = X) s =i °s = 0, fco is fixed up uniquely. 

The next theorem establishes that R(x, y, A) is indeed the kernel of the resolvent 
of L(X). 

Theorem 3.1. Let q{x) satisfy conditions (C.l) and (C.2) and let A^ be the 
simple zeroes of the minors irit(X). Then 

1. R ± (x,y, A) is an analytic function of X for X e C± having pole singularities 
at Xf € C±; 

2. R^(x, y, A) is a kernel of a bounded integral operator for ImX ^ 0; 

3. R(x, y, A) is uniformly bounded function for X G K and provides a kernel of 
an unbounded integral operator; 

4. R^(x,y, A) satisfy the equation: 

(3.6) L{X)R ± {x,y,X) = iS(x-y). 

Idea of the proof. 1. is obvious from the fact that x ± ( :E i^) are the 
FAS of L(X); 

2. Assume that ImA > and con side r the asymptotic behavior of R + (x,y, A) 
for x,y — * oo. From equations ( |2.9| ) we find that 

n 

(3.7) R+(x,y, A) = ^^x, X)e~ iXa ^-^Q+ p {x - y)£+(y, A) 

p=i 

Due to the fact that x + i x i ^) has triangular asymptotics for x — > oo and 
A G C+ and for the correct choice of Q + (x — y) (3.4) we check that the right 
hand side of (3.7) falls off exponentially for x — > oo and arbitrary choice of 
y. All other possibilities are treated analogously. 
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(3.8) 



For Agl the arguments of 2) can not be applied because the exponentials 
in the right hand side of (3.7) Im A = only oscillate. Thus we conclude that 
R ± (x, y, A) for A G M is only a bounded function and thus the corresponding 
operator R(X) is an unbounded integral operator. 

The proof of eq. (3.6) follows from the fact that L(X)x + (x, A) = and 
d@ ± (x - y) 



dx 



T&S(x-y). 



□ 



PROPOSITION 3.2. Let q{x) satisfy the co nditio ns (C.l) and (C.2), let 3 = 
{\f, j = 1, ■ ■ ■ , N} be the set of simple zeroes ( 2.2L ) of the minors mf(X) and let 
Ij < fco < Fj for all j — 1, . . . , N. Then the kernel of the resolvent R + (x, y, A) 
(resp. R~(x, y, A) ) has simple poles for X = X^ (resp. for A = Xj) with residues 
given by: 

(3.9a) Res^R ± (x,y,X) = \nf (x)) (mf (y)\ , 



(3.9b) 
(3.9c) 



n 



(*)) 



p j( x ))xtj( x ) U o\n ,j), 

Ujix)) 



K(y)| 



' 3 ^ " (™j 0*01% 0*0)' 

where A^ = pbj ± iUj and \q )( x ) 



( m j (y)| 



X^(x,Xf) 



i m j(y)\ n j(y)) 

(moj\lhjfiAv){i-Pi(v)), 



are the FAS corresponding to the 
potential qo satisfying (C.l) and (C.2) and whose set of simple zeroes is 3o = 
3\{A+A7}. 

Proof. Let Xo ( x , A) be the FAS of L (A) with potential qo(x); then Xq{x, A) 
are regular for A = A^. Now we apply the dressing method choosing A~^ as the 
locations of the singularities and construct the projector Pj(x) using the constant 
vectors |«oj) and (moj\. The normalizing factor uJ_(A) in the right hand side 
of (2.54) is a diagonal matrix that commutes with IIo. Then we insert X i x i A) = 
Uj(x,X)xo(x,X) in (3A) and note that the pole singularity of R + (x. y, X) at A = A+ 
(resp. R~(x,y,X) at A = Xj) can come up only from the factor uJ 1 (y : X) (resp. 
u(x,X)). To derive the ex press ions i n (|3.9|) one needs the explicit form of the 
projectors Pj(x) and Pj(y) ( 2.57 ) and ( |2.59| ). 

The right hand sides of (|3.9| ) do not vanish if the following conditions 

n oKj) ^ \n ,j), or n |n ,j) ^ 0, 

(3.10) 



(m ,j|n ^ (to , 



or 
or 



(TO 0j |n ^o. 



hold. In other words if ( 3.10 ) hold then the residues (|3.9| ) do not vanish, R ± {x, y, A) 
have simple poles at A = A^ and by definition A^ are discrete eigenvalues of L(X). 
Eq. ( |3 .10 ) is equivalent to the condition Ij < fco < Pj- Indeed violating this 
condition we get either (1 — n )|noj) = or n |noj) = and as a result - 
vanishing right hand sides in ( |3.9| ). 

To finish the proof one must check that from the definitions (3.9b) the relations 
( 2.68) ) follow. Besides |n^) and (m*| satisfy: 



(3.11) 



dx 



+ (q(x)-Xfj)\nf)=Q, 



d(mf 

dx 



-(mf\(q(x)-Xf) = 0, 
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where q(x) is given by ( |2.63 ). 



□ 



Corollary 3.3. The discrete spectrum of the Lax operator ( p.q j consists of 
the zeroes of the principal minors rnj'(X) for X G C+ and mj (A) for A G C_ 
provided the conditions (3.11) are satisfied. 

Now we can derive the completeness relation for the eigenfunctions of the Lax 
operator ( |l.5| ) by applying the contour integration method (see e.g. |27], Q) to 
the integral: 

(3.12) 2{x,y) = ^-.i dXR+(x,y,X)- -L <f dXR-{x,y,X), 

where the contours j± are shown on the Figure [l]. Skipping the details we get: 



n , , pea ( fco 

*(*-,)£-^ = ^/ ja E 

s=l s J -°° (s=l 



|x W+ (x,A))(x W+ (2/,A)| 



(3.13) 



■ £ lx [sI "(^A))(x w -(y,A)| 

AT 

2j £ k» { I n t (*)) ( m t (y)\-\ n 7 (*)) ( m 7 (y) I } 



3=1 

This relation ( 3.13] ) allows one to expand any vector- function (2(1)) G C" over 
the eigenfunctions of the system (1.5). Indeed, let us multiply (3.13) on the right 
by J I z(y)) and integrate over y. This gives: 

I *(*)> = £ fZo \x ls]+ (x, A)) ■ C(A) - E' s U, + i X W "(x, A)) • C(A)} 

(3-14) + Ef=i ^ (|4W> C+ - |n7(x)) C) , 
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where the expansion coefficients are of the form: 
(3.15) 

i>QO 

C±(A) = -i/ dx(x ls]± (x,X)\J\z(x)}, r ± (A) = -i/ dx(mf\j\z(x)). 



Remar k 3.4 . If q(x) ~ then x + ( x , ^) — X ( x , A) ^ exp(-iXJx) the set 3 is 
empty and ( 3.14 ) goes into the usual Fourier transform for the space C n . 

Remark 3.5. Here we used also the fact that all eigenvalues of J are non- 
vanishing. In the case when one (or several) of them vanishes we can prove com- 
pleteness of the eigenfunctions only in a certain subspace of C™. 

The resolvent for the CBCs is defined quite analogously: 
R(x,y, A) = R v (or, y, A), A e f2„, 
R u {x, y, A) = ixu{x, X)Q u {x - y)x„(x, A), 
(3.16) G»( Z ) = 0(-z)I%-9(z)(t-IlZ), ipj = £ £ ss , 

where Xu{ x iX) — £,u( x , X)e tXJx and k 0u is the number of positive eigenvalues of 
Im (A J) in the i^-th ordering. 

The following theorem is a specific case of one in 



Theorem 3.6. Let q{x) satisfy the conditions (C.l) and (C.2) and let 3 = 
U" = i (32p-i U 32p) where 

3ap-i = {A;^ 1 G fi2p_i, j = 1, . • . ,N} , 

(3.17) 3 2p = {\Ju p G fi 2p , j = l,...,iV}, 

are i/ie sets o/ zeroes and poles of the minors m v ^{\) in the sectors VL U . Then 

1. R u (x,y, A) is an analytic function of X for A £ 57„ having pole singularities 

at 3z,; 

2. i?„(x, y, A) is a kernel of a bounded integral operator for A € fij,; 

3. For A G /j/ U R u (x,y, A) is an uniformly bounded function which is a 
kernel of an unbounded integral operator; 

4. R u {x,y 1 A) satisfies the equation: 

(3.18) L(\)R v (x,y,\) = !S(x-y). 

The next natural step is to establish the structure of the singularities of R u (x, y , A) 
at the points of 3- This is done quite analogously by using the dressing factor ( 2.72 ). 
Note that in these matters the symmetry complicates the calculations. 

One of the effects of the Z„-symmetry is that the sets 3i/ are determined 
uniquely by 3i and 3o : 

(3.19) 3i = {Ate fi ls j = l,...,N}, 3o = {AJe02„, j = l,...,N}. 
The residue of R u (x, y, A) at the point A = X" can be cast into the form: 

Res Rx(x,y,X) = -2iImA+|n+(a;))(m+(x)|, 

A=A+ 
j 

(3.20) Res_R 2n (x,y,X) = 2iImXj\nj(x)){mj(x)\ : 
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where |n~~(x)) and (m^(a;)| are properly normalized eigenvectors of the Lax oper- 



ator corresponding to the eigenvalues A* G f2±i. The residues in the other sectors 



Q, v with v ^ 0, 1( mod 2n) are evaluated fro m (p.20| ) by employing eq. ( |2.31[) . 
Here we also have the analog of the condition ( |3.10 ). 

The derivation of the completeness relation of the eigenfunctions for CBCs 
with Z„-rcduction follows the same lines but needs some modifications. Instead of 
Z{x,y) ( 3.12| ) we should consider 

(3.21) 3(x,y)=J2 K ^, j dXR v (x,y,X), 
where the contours 7^ are defined by: 

(3.22) 72w _i = l a „_i U 7^_! U l 2u , l2u = hu U ITu+i U 72^ 



Here l v are the rays ( 2.24 ) oriented from to 00; 7^° is the 'infinite' arc Roe llfo 
with i?o > 1 and ir(i> — V)/n <ipo < irv/n; by overbar we denote the same contour 
with opposite orientation. Thus all the contours 721/-1 (resp. 72^) are positively 
(resp. negatively) oriented. 

Now we apply again the contour integration method and get two answers for 
Z(x,y). The first, according to Cauchy residue theorem is: 



2n N 



(3.23) Z(x,y) = ) ) Res R 2p+ i(x,y, X) + Res R 2p (x,y,X) \. 

P=l ] = 1 \ 3 3 / 

Integration along the contours taking into account that lim^oo x v ( x , A) = 1 gives: 
2,1 < r 



(3.24) Z{x,y) v ; . / di(J^(a;,tf ) A)-J2 I/ _i(a; ) i/,A)) + J-M(a;-i/). 

i/=i 717 

The completeness relation follows after equating both expressions and taking into 
account that (compare with ( |3.20| ) and ( 2.3l| )): 

Res R 2p+1 (x,y,X) = -2iImXf\nf p+1 \x)}(mf p+1) (x)\, 

j 

(3.25) Res R 2p (x,y, A) = 2ilm\+\nf p \x)){mf p) (x)\, 

A— ljp 

where \n^ p \x)), |nj 2p+1 - ) (x)) (resp. {ni^ p \x)\ , (mj 2p+1 '(a:)|) are properly nor- 
malized discrete eigenfunctions of the CBCs (1.5) (resp. of the adjoint CBCs ( |2.2| )) 
corresponding to the discrete eigenvalues Xjuj 2p and X^uo 2p . For the lack of space 
we can not provide all the details of the calculations. The final result is similar to 
the one for GZSs. Namely, any vector-function \z(x)) S C™ can be expanded over 
the eigenfunctions of the CBCs as follows: 

(3.26) \z{x)) = 

/ dX \ £ CA x )\x uAs] (x,x))- £ G,.W\x v ~ 1,M ( x > x )) 

v=l l " s < ko,„ s > ko.v 



N 2n 



LL toA l [£ d Mz) v ' + )-Cj\M*r'-)] > 

j=l u=l 
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where the expansion coefficients are given by: 

C+ (A) = 



(3.27) 



C S (A) 



dx{r> [s Kx,\)\j\z{x)), 

OO 

dx(r- Us] {x,X)\J\z(x)}, 



C, = -i 



dx{m v Ax)\ J\z(x)), 



C~- = 



dx{ra^- 1 {x)\J\z{x)). 



The completeness relations derived for GZSs and CBCs above can be viewed 
as the spectral decompositions for the generically non-self-adjoint operators L(X). 

Remark 3.7. The special case of a CBCs with Z„-symmetry is equivalent to 
n-th order scalar differential operator [11 1. Indeed, one can easily check that the 
system L ([j~5|), fll.SP can be written down as: 



(3.28) 



dx 



fe=i 



fc=i 



X(x,\)=0. 



After similarity transformation with uq = ■_ 1 u) s ° E a j goes into 



L X = ~u lLu oX 



d 
dx 



s=l 



s (x)E ss -\J2 E - 



s,s+l 



(3.29) 



4> s {x) 



k=l 



X(x,\) =0, 
A = iXcQui^ 1 / 2 , 



and can be rewritten as the scalar operator 

(3.30) L^xi = d n d n -i ■ ■ ■ d 2 d lX i(x, A)) = ~X n xi(x, A)), 

where dkX(x, A) = dX/dx + 4>k{x)X{x, A). If 4>k(x) are real functions (additional 
Z2-reduction of the type (2.6) ensures this) then L^ n ' is a self-adjoint operator. 

Remark 3.8. The author is aware that these type of derivations need addi- 
tional arguments to be made rigorous. One of the real difficulties is to find explicit 
conditions on the potential q(x) that are equivalent to the condition (C.2) or equiv- 
alently, to the conditions that m^(A) have only finite number of simple zeroes. 
Nevertheless there are situations (e.g., the reflectionless potentials) when all these 
conditions are fulfilled and all eigenfunctions of L(X) can be explicitly calculated. 
Another advantage of this approach is the possibility to apply it to Lax operators 
with more general dependence on A, e.g., quadratic or polynomial in A. 

The 'diagonal' of the resolvent 

By the diagonal of t he re solvent one usually means R(x, y, A) evaluated at y = x. 
However the definition (3^3) is not continuous for y = x and needs regularization. 
The simplest possibility is to consider as the diagonal of the resolvent: 

R{x, X) = ^ (R(x + 0, x, A) + R(x, x + Q, A)) . 

In fact we will consider as the a somewhat more general expression: 

(3.31) R P (x,X)=ix ± (x,X)Px ± (x,X), 
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where P is a constant diagonal matrix. Obviously Rp(x, A) satisfies 
,dRp(x, A) 



(3.32) 



dx 



+ [q(x) - AJ, R P {x, A)] ee [L(X), R P (x, A)] = 0. 



Thus Rp(x, A) belongs to the kernel of the operator [L(X), •]. Due to the fact that 
X (x, A) is the FAS and satisfies a RHP with canonical normalization we find: 



(3.33) 



R P (x,X) =iP + J2R ( p\x)X- k . 



k=l 



The coefficients R^ (x) can be expressed through q(x) using the recursion relations 
generalizing the ones of AKNS 18, 37, 27|. These relations are solved by the 
recursion operators A± which have the form: 

dX 



(3.34) A±A = ad7 J 



dx 



P ([q x ),X(x)]+i 



q[x), I dy[q{y),X(y)} 



where Pq is the proje ctor onto the off-diagonal part of the matrix PqX = Xq, the 
matrix X(x) in ( 3.34 ) satisfies X = P X and 



(ad - l X { ) 



a,i — ou 



The coefficients R^ (x) can be expressed in compact form through A± as follows: 
(3.35) R k P +li = A±R P i = -k k ± &&- J 1 [iP 1 q{x 1 t)] 1 



(3.36) 



4 fc)d = i 



dy(l-P )([q(y,t),R$ )f \ 



lira R 

x— *±oo 



(fe)d, 



x,t). 



Quite naturally these coefficients, or rather the diagonal of the resolvent gen- 
erates 0, 0, |J|: 

- the class of NLEE. Given the dispersion law, e.g., /(A) = X N P of the NLEE 
we can write down the equation itself by: 



(3.37) 



dq . 



,p , 



dx 



+ P ({q(x,t),R { p N \x,t)}=0. 



- the corresponding Lax representations, or in other words, the M-operators 
for each of these NLEE as follows: 



(3.38) 



V^(x,X) = ^R^(x)X 



x N-k 



k=0 



the integrals of motion of the corresponding NLEE. This follows from 



Theorem 3.9 (|18[). The quantities 
(3.39) R^ w (x,X)=i X ± (x,X)IL^x ± (x,X), II« = V E ss - -1, 
satisfy the relations 



s=l 



(3.40) 



/ dxtr (r£ w (x, A) J - ffl (fc) J 

•J — oo ^ 



dX 



where S)^(A) is defined by (2.4t )■ 
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Combined with the (3.34) we can deduce that the diagonal of the resolvent and 
the recursion operator 



(3.41) 



(A± - X)R$,(x, A) = i[P, ad ^q{x) 



directly reproduce the generating functionals of the conserved quantities. 

The termin 'squared' solutions and recursion operator do not reflect properly 
the algebraic properties of these objects. The recursion operators A± can be under- 
stood as the Lax operator L(X) in the adjoint representation. One of the definitions 
of the adjoint representation means that we s hould replace each element U(x, A) € g 
by ad jjr Xj \y — [U(x, A), •]. Therefore due to ( |3 .32 ) we can view the diagonal of the 
resolvent Rp(x,X) as the eigenfunction of L(X) in the adjoint representation. It 
remains to project out the kernel of ad j in order to derive A± from L(X). 

The 'squared' solutions are eigenfunctions of A± and belong to a linear space, 
which is the co-adjoint orbit of q* + determined by J. The gauge covariant way to 
introduce them involves the FAS of L(X) and is: 

(3.42) 

e%{x,X) = P ( X ± (x,X)E ij x ± (x,X)) , hf(x,X) = P ( X ± (x,X)H jX ± (x,X)) 

where % (a;, A) are the FAS of L(X) GZSs. The similarity transformation by 
X {x, A) is the adjoint action of the group 25 on the algebra g; therefore e^(x, A) 
and hj (x, A) are elements again of q. The projection IIo = ad7 1 adj is a natural 
linear operator on g. Besides the 'squared' solutions are analytic functions of A 
having both poles and zeroes at A^ . 

More detailed analysis based on the Wronskian relations reveals several other 
important aspects [36, 19, 30 of the 'squared' solutions of GZSs. First, the sets 



{e+{x, A), e^O, A)}, e±. k (x), e jiik (x), e±. k (x), e^x), 
and 

{e±(x, X),er j ( x , A)}, e± k (x),e^. k {x), e±. k (x),eZ k (x), 



i<j,k = l,...N} 



i<j,k = l,...N} 



form complete sets of functions on SOT that realize the mapping 9Jt <-> T. Here by 
ef tk (x) and e^. fe (x) we have denoted: 



-ij-k 



(x) 



defjx,X) 



dX 



x=xf 



Second, it is possible to expand the potential [P,a,dj 1 q(x,t)] and its varia- 
tion a,d^ r 1 Sq(x) over each of the complete sets shown above. The corresponding 
expansion coefficients are expressed through T and their variations. These facts 
constitute the grounds on which one can show that the ISM can be understood as 
a generalized Fourier transform. The important difference as compare to the stan- 
dard Fourier transform is in the fact that the operator L (as well as the operators 
A±) allows for discrete eigenvalues. Therefore the completeness relations involve 
both integrals along the continuous spectrum and sum over the discrete eigenvalues. 
In the usual Fourier transform the discrete eigenvalues are absent. 
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Hamiltonian properties of the NLEE 

Here we briefly formulate the Hamiltonian properties of the NLEE paying more 
attention to its algebraic structure. This has been widely studied problem, see 
(|, Hi], |l4], 0, [uj[ [l3|, fl|, HH, and the numerous references therein. 

In doing so we follow mainly the ideas of [39] with a natural generalization 
from sl(2) to sZ(n)-algebras. The main idea in these papers is the possibility to 
write down the Lax equation (1.4) in explicitly Hamiltonian form as the co-adjoint 
action of g on its dual g* . Obviously any non-trivial grading in g (resp. g, g) will 
reflect into a corresponding grading of the dual algebra g* (resp. g*, §*). 

Below we will need also the Cartan-Weyl basis of sl(n). Choosing for defmite- 
ness the typical n x n representation we fix it up by: 



(3.43) fj = I.e. {H, = E u - E, 



i+l,i+l> 



I = 1, 



1} 



{Eij, i j] 



As invariant bilinear form we can use (X,Y) = trpTY"). Then the commutation 
relations can be written in the form: 



(3.44) 



[Hi, Ejk] 
[Ejk,Eki] 



(e, - 
Ejt, 

k-l 



e-i+l,Cj - e k )E jk , j ^ k, 
[E jk ,E lj \=-E lk , l^j, 



By e k above we mean an orthonormal basis in the n-dimensional Euclidean space 



with a standard scalar product: (ej,e k ) 



Ojk- 



Those, who are familiar with Lie 



algebras will recognize ej — e^+i as the simple roots of sl(n) and the set of ej — e k , 
j ' ^ k as the root system of sl(n). 

If C = 1 (i.e. with the trivial grading) each of the matrices U k (x) in ( 1.19] ) is 
of generic form: 



(3.45) 



U k (x) 



n-l 
i=i 



E 



(fc) P 



The coefficients Uj k \x), u^(x) can be viewed as linear functionals on u k (x) and 
thus they belong to g*. Using the bilinear form (1.27) they can be interpreted as 
linear functionals on g and thus as elements also of g* . The algebraic structure on 
g* can be introduced in analogy with the commutation relations (3.44), namely: 

m} 

> v 

j P 



, \ (s+m—p) 

{ei-e i+1 ,ej - e k )u) J+k 



(x)8(x - y), 



(3.46) |i 



(«) 



,( m ) 



'i,i+k( X )' U i+k,j * 



(s+m—p) 



(x)6(x - y), 



i+k-l 

u ( s + m ~p) 

l=i 



n 



(x)S(x -y) + iS s+m , P S'(x - y). 



The derivation of these relations follows |39j in a rather straightforward manner; 
though a bit tedious, it can be generalized also to any simple Lie algebra. 

Note that if p = — 1 then the term with 8'{x — y) disappears and the Poisson 
brackets (3.46) become ultralocal. Then we can rewrite them in a compact form 
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using the classical r-matrix [13]: 



(3.47) \ U(x, A) $ U(y, fi) } = [r(X - fi), U(x, A) ® 1 + 1 <g> U(x, (j,)]S(x - y), 



(3.48) 



r(A - ,i) = 



n„ 



A 



The left hand side of (|3.47 ) has the structure of the usual tensor product oftixn 
matrices, but instead of taking the product one should rather take the Poisson 
bracket between the corresponding matrix elements of U(x, A) and U(y, //). 

The relations (3.47) are local in the sense that for the evaluation of the left 
hand side of (3.47) we need to use only the Poisson brackets between the matrix 
elements of U(x,X) and do not need the boundary conditions on the potentials. 
The effectiveness of the r-matrix, when it exists, is in the possibility to evaluate 
the Poisson brackets between the matrix elements of the scattering matrix T(A). 
To do this we need to 'integrate' (3.47) which needs to take into account also the 
boundary conditions. For periodic boundary conditions on q{x) this gives: 



T(A)®T(/M = [r(A- M ),T(A)®2»]. 



(3.49) 



For vanishing boundary conditions on q(x) and J = J* the calculations need some 
additional considerations with the result (see [13]): 



T(A)®T(/i) \ = r + (\ - fi)T{\) ®T(fj,) - T(\) ®T(fi)r-{\ - fx), 



(3.50) 



r± (A - fi) = ^2 Ejj ® Ejj T inS(\ — fi) ^ J5y (g> E 3 



From both relations ( [3.4£ ) and ( 3.5C ) there follo ws th at the principal minors m^(A) 
commute with respect to the Poisson brackets ( fc.46 ) [ 19 1 



i.e. 



(3.51) 



{S) fc (A),3) J -(M)}-i=a 



Since £>fe(A) are the generating functionals of integrals of motion S^' (see eq. 



(2.51)), then eq. (3-51) means that all these integrals are in involution with respect 
to these Poisson brackets. 

The Z„-symmetry may modify substantially some of the above results. Indeed, 
it can be viewed as a set of constraints on the phase space 9Jt and on the generic 
Poisson brackets (3.46). Then one should evaluate the corresponding Dirac brackets 
on the reduced phase space. However in the case of the Z„-NLS equation ( [1.3] ) with 
Lax operator L given by (1.5), (1.8) somewhat surprisingly the approach of ]39| 
gives us directly the correct answer. If we define ipj(x,t) as linear functionals of 
U{x, i, A) = q(x, t) — A J by: 



(3.52) fa (x, t) = -tr (U(x, t, \)K n ~ j ) , 

and make use of fll.Sj ) then the set of Poisson brackets in ( |3.46 ) simplify to 

(3.53) {ipj(x,t),ip k (y,t)} =S k+j - n S'(x-y). 
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Together with the Hamiltonian H — lo 2 m[ 2 1 ( 2.53 ) they provide the Hamiltonian 
formulation of (1.3). Unfortunately this Poisson brackets are not ultra- local and 



the corresponding Lax operator does not allow classical r-matrix of the form (3.47) 
For the affine Toda chain (1.2) the simplest Poisson brackets are provided by: 

dx 



(3.54) 



kiV,t) > = 5 kj 5{x - y). 



The correspondin g Lax operator (|L!^) unlike the previous case allows classical r- 
matrix satisfying ( [3.47j ) which however has more complicated dependence on A — /i; 
it is known as the trigonometric r-matrix [38]. 

Another special property of the Z n -symmetric CBCs concerns the existence 
of the so-called symplectic basis The el emen ts of these bases are special 

linear combinations of the 'squared solutions' ( |3.42 ) which are also complete in 
dJl and which are such that the expansion coefficients of 5q(x,t) over it produce 
the variations of the action-angle variables of the corresponding set of NLEE. In 
|25| this basis was worked out for the Zakharov-Shabat system related to the sl(2) 
algebra. For GZSs related to algebras of higher rank such basis is yet unknown 
although it must exist since the action-angle variables for them are known [^, Q . 

For the Z n -symmetric CBCs the construction of the symplectic basis is very 



much like the one in 25 1 due to the fact that the subalgebras q u related to each of 



the rays l v are direct sums of sl(2) subalgebras. It is a complete set of functions on 
the phase space of the corresponding Z„-symmetric NLEE (1.1) and (p"^). Skip- 
ping the details we just give the explicit expressions for the set 21 of action-angle 
variables of the Z„-NLS equation in terms of the scattering data of its Lax operator. 
Obviously 21 will consists of two sets of functions 21 = 2lo U 2li each set defined on 
the ray Iq and l\ respectively: 



2l = {Tr lJ (X),K lJ (X), 

2tl = {TTij(X),Kij(X), 



A £ l , 
A G h, 



i + j 



2( 
1( 



where 



(3.55) 



r «(A) = In I 

7T 



pti pi- 



rn 



'2 bZ (A)' 



mod n)} , 
mod n)} , 

piM = 



KM) 

4(A)' 



bfj(X) were introduced in ( 2.30| ). 



and the coefficients a^(A) 

Quite analogous are the expressions for the action-angle variables for the two- 
dimensional affine Toda chain provided we use the scattering data of the Lax op- 
erator fit. i op. 



4. Conclusion 

The restricted space did not allow us to give more details or explanations on 
these and related problems. We only mention some of them below. 

One such important to our mind result is the interpretation of the ISM as a 
generalized Fourier transform. In its derivation for the GZSs and CBCs [27, [l9| , pTo| 
both algebraic methods and analytic ones were used. As a result the pair-wise 
equivalence of the symplectic structures in the hierarchy becomes obvious. 

The approach based on the Kac-Moody algebras is a natural basis for the 
Hamiltonian hierarchies. If one can derive a bi-Hamiltonian formulation of a given 



:!() 



VLADIMIR S. GERDJIKOV 



NLEE then there is a whole hierarchy of them related by a recursion operator A. 



Here we mention the paper [15 where the operator A was derived as the 'ratio' of 
two such Hamiltonian structures for the N-w&ve equations. The result, of course 
coincides with the natural expression for A obtained with the AKNS recursion 



method and whose spectral theory was constructed by other means in [27, 18 



The method based on the diagonal of the resolvent of the Lax operator started 



by Gel'fand and Dickey |17, 10 1 can be viewed also as a formal algebraic one. 
The authors studied by algebraic means the ring of operators, commuting with L. 
They expressed most of the quantities, including the diagonal of the resolvent of 
L, as series over fractional powers of L and did not investigate the existence and 



convergence of these series. Once identified with the expression (3.31) in terms of 
the FAS these problems find their natural and positive solution. 

Besides the classical r-matrix corresponding to the ultralocal Poisson brackets 
there exist also dynamical r-matrices depending on the fields qij(x) in the NLEE. 
One of the problems, that is still not solved is to find the interrelation between the 
dynamical r-matrices, r and the recursion operator A. 

Finally, we should mention that both approaches have been further generalized. 
For example, the analytic approach was generalized from a local RHP to a nonlocal 



RHP and to <9-bar problem (also local and nonlocal), see Jy, 50, 37 1. This allowed 
to treat NLEE of soliton type in 2 + 1 dimensions. 

Another direction is to study Lax operators with more general A-dependencc 



such as polynomial, or rational [51 



Obviously all results concerning spectral decompositions can be formulated in a 
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gauge covariant way thus allowing to treat also gauge equivalent NLEE [28, 29, 

The algebraic approach was also generalized to use as a basis infinite dimen- 
sional algebras such as Virasoro algebra, W\+ao etc. which lead to the important 
construction of the Japanese r-function and its relation to the soliton theory, see 



32, 14 



Thus we just outlined the beginning of all this and so it is time to stop. 
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Appendix A. Gauss decompositions 

The Gauss decompositions mentioned above have natural group-theoretical in- 
terpretation and can be generalized to any semi-simple Lie algebra. It is well 
known that if given group element allows Gauss decompositions then its factors are 
uniquely determined. Below we write down the explicit expressions for the matrix 
elements of T ± (A), S ± (X), D ± (\) through the matrix elements of T(A): 

(A.l) T-.(A) = 

(A.2) fr(A) = 



in 



1, 2, p 

+ (A) 11-2 ., \.j 



T(A) 



i-iy+p fi, 2 „ j \w 



i(A) \l, 2, p, j-1 J 



T(A) 
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(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 
where 

(A.9) 



55(A) = 
S+(A) = 



s-(x) 



1, 2, 



rot (A) I 1, 2, 



1, 2, 



n-j+l 
(_l)P+j 



P, 



P, 



i - 1, i 1 

3 -hP) 



, j-1 

, 3 



CO 

T(A) 



1, 



, +1 (A) U J + 



J, j + 1 
■ (X) {j + 1,3 + 2 



(-l) p+j 



j + J + 2 

i, j + 1 
i, j + 



+i(A) IP, J + 



T(A) 

rc. — 1, n 
n — 1, n 

■ j P> • • • , 
■ , P, • • • , 

■ , Pj • • • , 
. , p, 

71—1, n 
n — 1, n 



(n-j+l) 

T(A) 
(n-j) 



T(A) 

T(A) 
(n-j+l) 



T(A) 



«1, «2, 
Jl) 5*2) 



3k 



(fe) 

T(A) 



det 



T 

J- V 



T- ■ T 



T- 
T 

■ J- 7.' 



T- T- ■ 



ikjk 

i k and the columns 



is the minor of order k of T(A) formed by the rows i\, £2, ■ 
ji, J2, ... , jk\ by p we mean that p is missing. 

From the formulae above we arrive to the following 

COROLLARY A.l. In order that the group element T(A) £ SL(n,C) allows the 
first (resp. the second) Gauss decomposition ( 2.1(\ ) is necessary and sufficient that 
all upper- (resp. lower-) principle minors mi" (A) (resp. mZ{X)) are not vanishing. 

These formulae hold true also if we need to construct the Gauss decomposition 
of an element of the orthogonal SO(n) group. Here we just note that if T(A) G 
SO(n) then 



(A.10) 
where 

(A.ll) 



5 (T(A)) T 5 - 1 =T- 1 (A), 



So = Y^(-l) k+1 {E k>n+1 _ k + E n+1 - k . k ), if n = 2n , 
fe=i 



S = ^(-l) fe+1 (£fc,«+i-fc + E n +i-k,k) + {-l) n °E no+hno+1 , if n = 2n + 1. 



k=l 



One can check that if T (A) sa tisfies ( |A.10|) then each of the factors T ± (A), 5 ,± (A) 
and D ± (A) also satisfy ( A. 10 ) and thus belong to the same group (3. In addition 
we have the following interrelations between the principal minors of T(A): 



(A.12) 



mf(X)=m±_ j (X), 
m±(A)=m± _,(A), 



for SO(n), 
for 5P(n), 
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Appendix B. Dispersion relations for Sfc(A) and \nm^ k (X) 
Let us introduce the functions f^{X): 

' 'A -At 



#(A) = 



Rk(X) ' 



N 



f-_ k (X) = R k (X)m-_ k (X), 



= n 



a- a: 



which like m^(A) are: i) analytic for A e C±; ii) satisfy lim^oo fjtW = 1- Besides, 
/if (A) have no zeroes in their regions of analyticity and therefore the functions 
ln/k(A) are analytic for A G C± and tend to for A — > oo. This allows one to apply 
the Plemelji-Sokhotzky formula with the result: 



(B.l) 

where 
(B.2) 



Sfc(A) 



1 

2vri 



dfi 
jJL — A 



ln/"_ fc (A) 



A e C+ 
A g C_. 



It remains to insert the above definitions of /r(A) into (B.l) and (B.2) to derive 
Eqs. (fEIsj), (|2~46|). 



The dispersion relation (2.5C) is derived analogously treating the integral 



(B.3) 3(A) = J2 



2?ri 



dfi 
_ /i-A 



In < m. 



n N 

»nn 



/x — A 



/i — A 



with Ae!!„ and the contours 7, as in (3.22). 
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